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THE /-FUNCTION DISTRIBUTION AND ITS EXTENSIONS 


P. VELLAISAMY AND K. K. KATARIA 

Abstract. In this paper we introduce a new probability distribution on (0, oo), associated 
with the /-function, namely, the /-function distribution. This distribution generalizes 
several known distributions with positive support. It is also shown that the distribution of 
products, quotients and powers of independent /-function variates are /-function variates. 
Another distribution called the /-function inverse Gaussian distribution is also introduced 
and studied. 


1. Introduction 


Springer and Thompson jUj showed that the probability density functions (pdf’s) of 
products and quotients of independent beta, gamma and central Gaussian random variables 
(rv’s) can be expressed in terms of the Meijer G -function (see [HI], p. 21). Matliai and 
Saxena [T^j discussed the distribution of the product of two independent rv’s whose pdf 
can be expressed as the product of two //-functions; see [g], [l5[ for definition and more 
details on the //-function. Carter and Springer Q introduced the //-function distribution 
which includes, among others, the gamma, beta, Weilbull, chi-square, exponential and the 
half-normal distribution as particular cases. Inayat-Hussain f§] generalized the //-function, 
namely, to the //-function. The generalized Riemann zeta function, the polylogarithmic 
function of complex order and the exact partition function of the Gaussian free energy 
model in statistical mechanics are some special cases of the //-function, which are not 
particular cases of the //-function. Rathie [l8[ introduced the /-function, which includes 
the //-function as a special case. The /-function is represented by the following Mellin- 
Barnes type contour integral 


where 


!{*) = C" 


(flj, Aj, C%i)l } p 

C bj , Bj, /3y)i,g 
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X(s)z S ds, 
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riLn+l r ai (Oi + A i S ) n?=m+l ^ (! - bj - BjS) 


( 1 . 1 ) 


( 1 . 2 ) 


An equivalent definition can be obtained on substituting w = — s in (II. 1)1 . In the above 
definition, z ^ 0 ,m,n,p,q are integers satisfying 0 < m < q and 0 < n < p with 
cui,Ai > 0 for i = 1,2 and (3j. B 3 > 0 for j = 1,2 Also, af s and bj's are 

complex numbers such that no singularity of (bj + Bjs) coincides with any singularity 

of T Qi (1 — a* — Ais). An empty product is to be interpreted as unity. The path of integra¬ 
tion C, in the complex s-plane runs from c — ioo to c + ioo for some real number c such 
that the singularity of Y 8j (bj + Bjs) lie entirely to the left of the path and the singularity 
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of T Qi (1 — a,i — A t .s) lie entirely to the right of the path. The conditions for convergence of 
the integral involved in (II.ip are as follows (see |l8j). Let 


q p 

V = PjBj — 'y ^ QijAi, 

j=l i =1 

n = Y a i~Y Pi’ 

i= 1 ' ' j =1 ' ' 

m q n p 

a = yPj b j - Y Pj B j + Y aiAi ~ ^2 

j= 1 j=m+ 1 2—1 2=72+1 


where m, n, a,i, +, a* and bj, Bj,/3j appear in the definition of the /-function. The /-function 
is analytic if V > 0 and the integral in (II.ip converges absolutely if | arg(z)| < A|, where 
A > 0. Also, if | arg(z)| = A^ with A > 0, then it converges absolutely under the following 
conditions: (i) V = 0 and £7 < —1, (ii) |V| ^ 0 with s = cr + it, where a, t 6 M and are 
such that for \t\ —> oo we have 17 + aV < —1. The //-function follows as a particular case 
when a* = 1 for i — n + 1,... ,p and + = 1 for j = 1, 2,. .. ; m in (II.ip . Recently, the 

and its applications 
, H3). The effect of 
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extension of the /-function to several complex variables (see [11[ 
to wireless communication were studied by several authors (see 
Marichev-Saigo-Maeda fractional operators on the /-function has been recently studied by 
Kataria and Vellaisamy 


see 


10 


The main purpose of this paper is to introduce a new probability distribution based on 
the /-function, which we call the /-function distribution. In Section 2, some preliminary 
results are stated that will be used in subsequent sections. The Mellin and Laplace trans¬ 
forms of the /-function are evaluated in Section 3. In Section 4, the /-function distribution 
is introduced and some results related to the distribution of products, quotients and pow¬ 
ers of independent /-function variates are derived. Also, the /-function inverse Gaussian 
distribution is introduced and studied in Section 5. 


2. Preliminaries 


The /-function is connected to the generalized hypergeo metric function p F q , the Meijer’s 
G-function, the generalized Wright function p i]) q , the Fox’s //-function and the //-function 
by the following relationships: 
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The following identities will be used later: 
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(a + 6 + 1,1,1) 
( 6 , 1 , 1 ) 


= z b (5-z) a , \z\ < 6. 


Also, from the definition of the /-function, the following properties are immediate. 
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( 2 . 1 ) 

( 2 . 2 ) 
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( 2 . 6 ) 


Note that whether a product of two /-functions is another /-function is still an open 
problem. 


3. The Mellin and Laplace transform of the /-Function 

In this section, we evaluate the Mellin and Laplace transform of the /-function. The 
Mellin and Laplace transform of a continuous positive rv X with pdf fx is defined by 


and 


M fx (s) = E[X s ~ l ] = / x s ~ 1 f x (x)dx (3.1) 

Jo 
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£ fx (r) = E[e- rX ] = / e~ rx f x (x)dx, (3.2) 

■Jo 

respectively. Under some suitable restrictions 22j on A4f x (s), there exists an inversion 
integral 

™ = 2 b 

The following results are due to Epstein 
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(iii) Let X±,X 2 be two positive continuous independent rv’s. The Mellin transform of the 
rv Y = Xi/X 2 is 

M, r (s)= M, x js)M lx fi-s). (3.6) 

(iv) Let X il X 2 ,..., X n be n positive continuous independent rv’s. Then the Mellin trans¬ 
form of the rv Y = YYj=i Xj is 

n 

M lr (s) = \[M, Xl (s). (3.7) 

3 = 1 


3.1. The Mellin transform. Assuming the convergence of the integral involved in defi¬ 
nition of the /-function, for a > 0 and <5 7 ^ 0 , we have 
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Substituting s = 
/-function (13.81) is 
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x(w)S w x aw dw. 
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aw in the above integral and using (13. dp . the Mellin transform of the 


= AAh. ( 3 . 9 ) 

The above result can be extended to obtain the Mellin transform of product of two I- 
functions. We state the following useful result in the form of a lemma. 


Lemma 3.1. The Mellin transform, A / lj(Sx cr )i(r}x^)(s), of the product of 


I{8x°) = / 


m,n 

PA 


Sx° 


(&Z5 ^z)l,p 

(Pji Bj , Pj)l,q 


and I(rjx' J ) = /*;* 


r)x H 


(ci, Q, 7i)i, u 
(dj, Dj, Pj)i t i 


for a > 0, p > 0 and 5,r) 7 ^ 0 is equal to 


1 7-n+fc,m+Z 

V 

(q? ? Tz) 1,Z 

(l - bq 

— a 1 BjS,pa 1 Bj,(3j 

)l,q 

(q, / y«)/+i,w 



(' dj ) Dj , Pj)i,k 

(1 - a t 

- a~~ 1 A i s, pa~ l A i} af 

l,p 

(dj , Dj , pj ) /c-|-l 5 t; 


Proof. From (13.11) . we have 
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We assume the conditions for absolute convergence of integral involved in above equation 
holds. By interchanging the order of integrals, we get 
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The proof is complete by using definition of the /-function (11.11) . 
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3.2. The Laplace transform. For simplicity, we use an equivalent definition of the I- 
function to determine its Laplace transform. Moreover, the conditions for absolute conver¬ 
gence of integral involved in the definition of I( K 5x a ) is assumed. From (13.21) . we have 

poo 

£i{ 5 x°)(r) = / e~ rx I(5x a )dx 


r»oo g—rx 


2tt i 


'c 


x(—s)(8x rT ) s ds dx 


4/+-^ I 


x as e~ rx dxds 


1 1 
r 27rz 


r(l + as)x(-s) ( —^ ds 


'c 


1 7-771,71+1 

~ J p+U 


r cr —1 

' 7-771,77+1 

£ 1 P+ 1,9 


n (J~ 1 

t7i+ 1,m 


( 0 , ( 7 , 1 ) 

(Pj 1 Bj, /3j)l,q 


(cr, (j, 1) (a* T Aj, ylj, c+gp 

{bj T Bj , -Bj, (3j ) i )9 


(1 Bj, Bj, /3j ) 1 ,<j 

(1-cr, a, 1) (1 - a* -+, Agt+gp 


(using ([276]) ) 

, (3.10) 


where the last step follows from (12.41) . 


4. The /-Function Distribution 

Here, we introduce a new probability distribution associated with the /-function. Let 
X > 0 be a rv with pdf 

( r 
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where <7 > 0, <5 7 ^ 0 and (a,, +, c+g p , ( bj , B 7 , /3,)g 9 , m, n must conform to those restrictions 
in the definition of the /-function (II. ip . such that fx( x ) > 0. Also, k 0 is such that 
r fx( x ) dx = l. 

(a) Moments. The r-th moment for the /-function distribution can be obtained using the 
Mcllin integral transform as follows. 

d'r = M fx (r + l) 
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using (13.9 ft and y(.) is defined in (11.2ft . The normalizing constant of the /-function distri¬ 
bution is obtained from // 0 = 1 as 


k 0 = 
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(b) The characteristic function. The characteristic function of the /-function distribution 
is given by 
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Hence, the moment generating function is 
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(c) The cumulative distribution function (cdf). The Mcllin transform of the cdf Fx(x) of 
a rv X is connected to its pdf fx (x) by the following relationship (see [20], p. 99). 
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The hazard rate function Ax(x) is given by 
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(d) Order statistics. Let Xp) < X( 2 ) < • • • < Xpv) denotes the order statistics of a random 
sample X\,X 2 ,... ,X^ of size N drawn from the /-function distribution. The pdf of the 
j'-tli order statistics is given by 

= (J3T) {FxWY' 1 {Fx(x )} N ~ J . 

Therefore, the pdf of the minimum and the maximum order statistics is given by 
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4.1. Some special cases. The density of the //-function distribution (see jij), 
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follows as special case of the /-function distribution (14.ip when a = 1, a^s and /3j , s for all 
i and j equals unity, ffence, the other standard non-negative distributions are particular 
cases. We next give the densities of the well-known distributions in terms of the /-function. 
Let X be a rv which follows 
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(iv) the Weibull distribution with parameters 6, A > 0. Then 

x > 0. (4.6) 

(v) the Rayleigh distribution (put 6 = 2 and A = 2u 2 in (14.61) 1 with parameter v > 0. Then 
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(viii) the half-Cauchy distribution with parameter A > 0. Then 
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(x) the beta distribution of second kind with parameters 9, A > 0. Then 
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(xi) the power function distribution with parameter 6 > 0. Then 
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(xv) the F distribution with parameters 9, A > 0. Then 
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The above relationships can easily be obtained using identities fl2.1j) - (12.3j) and properties 
of the /-function fl2.4jl - fl2.6p . 

4.2. Properties of the /-function distribution. Next we show that the class of /- 
function variates is closed under multiplication and quotients of independent /-function 
variates. Moreover, this class is also closed under positive scalar multiplication and rational 
power of the /-function variate. 

Theorem 4.1. Let Xl,X 2 , ..., Xjy be N independent /-function rv’s with pdf’s fx 1 (^ 1 ), 
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Now using the inverse Mcllin transform (13.3p . the pdf is given by 
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I Xi{°i 1 s) I v II S ° 


ds, 


1=1 


□ 


Example 4.1. Let Xi,X 2 ,..., X n be n independent rv’s such that X, follows half-normal 
distribution for % — 1 , 2 ,..., m and beta distribution for % — m + 1 ,..., n, with pdf’s 


»? 

2 o? 


V2 

fxi (a?t) = - 7= e 2at > x i> i = 1 > 2 , • • •, m, 


(Tie'll 


and fxX x i) = p f p 5 ! X T 1 (l~x i ) bi \ 0<x t <l, i = m + l,...,n, 

T(ai)T(bi) 

where the parameters ( 7 j,aj+ > 0. The corresponding Mellin transforms (13.1 ft are given 
by 


M fx (s) = 


vi<Ti (vii) r (s) > I - 

i = m + 1 ,..., n. 


r(cij)r(cii+6i—i+s) 5 


Let X = niU*, Using ( 13.71) . its Mellin transform is 

r(fli + frj)r(aj — i + s) 

2y/ r(a i )r(a i + b t - 1 + s) 


^w=r- n 


(4.2) 


Using (13.3p and (14.2p . the pdf of rv Y is 


friy) = 


k 
2n i 


r»c+zoo 


n 


i=m -\-1 


r (&i — 1 ~b s) 
r(ai + bi — 1 + s) ^ ^ 


n 


ds 


_ 7 jn— m+1,0 

n—m,n— ra+1 


o- — 

2 2 y 

nm 

i=l< 7 i 

_TT m l rr 71 r(ai+b^) 


(cti ~\~ 1 5 1} l)m+l,n 

(o, 2? 1 ? 1 ? l)m+l,n 


(4.3) 
y > o, 


where k = Jl™ i 75=7 TYi=m+i r \a ) ■ Note that the integral (14. 3 b can also be represented 
in terms of the Fox’s //-function as 


fy(y) = 


n n r (aj+bj) 

i=m+1 T(ai) rjnO 

~ - n„ 


(2 ^Yu=i°i 


2 2 y 

ifli b{ l)m+l,n 

TT m 

n,=i^ 

(0 5 2)1?^ 1) 1)771+1,71 


, y > 0 . 


Example 4.2. Let Ah, Ah be two independent rv’s such that Ah is uniformly distributed 
(I4.10|) on (0,1) and X 2 follows gamma distribution (14 . 4 p with shape parameter 9 = 2. Using 
Theorem 14.11 the pdf of rv Y = Ah Ah is 


fy(v) = +? 


y 


A 

"0 

1 _ 


, y > 0 , 


which is an exponential distribution with parameter A (see joj, p. 306). 
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Theorem 4.2. Let X be an /-function rv with pdf (14. ip . 
(a) The pdf of the rv Y = aX, a > 0, is given by 


fy(y) = 


6° 


rm,n 


1 

b 

*-© 

I_ 

(J -^4z5^z)l,p 

a 

(bj,a Bj,f3j)i !q 


ax (O 

(b) The pdf of the rv Y = X r for r ^ 0 rational is given by 


, y> o. 


fy(y) = 


X(a-i) ™ 
when r > 0 and 


5 ra ~ y 


(di-rcr 1 Ai + cr 1 A i ,ra 1 A i ,a i ) ltP 
(bj - ra^Bj + a^Bj, ra^Bj, /3j)i, q 


, V > 0 , 


fr(y) = 


5 r °- y 


X (a-i) ™ 
when r < 0. 

Proof, (a) From (13.41) . we have 


(1 - bj + ra l Bj - a 1 B j , -ra 1 B j ,/3 j ) hq 
(1 - a,i + m _1 + - a~ l A t , -rcr _1 +, af) hp 


y > c 


Mf Y (s) = a s 1 M fx {s) = 


a s l 5° 1 x{a V) 


Using (13.3j) . the pdf of Y is 


xa- 1 i r c+ioo / Sja- 1 ' 

fr(y) = — ( -i\o ■ / xO -1 s) -- ) ds , y > 0, 

ax((J Y 2m J c _ ioo 


which gives (14.41) on using (11.11) . 
(b) From (I3.5[) . we have 


M f Y (s) = M fx (rs - r + 1) = 


S ra 1 x( CT Vs ~ cr V + a *) 


Therefore, 


rc+ioo 


fy{y) = 


— 1 _ _ —1 i 1 \ / , era 


x( a V 2ni J c _t 


x(cr Vs — a r + a 


ds. 


When r > 0, the lhs of (14.71) on using (II. ip gives (14.51) . When r < 0, put t = —r 
for t > 0, 


fy{y) = 


5~ ta 1 1 


rc+ioo 


X(a V 2m J c _ ioo 

which on using (11.11) gives (14.61) . 


x(—cr Vs + a V + a x ) (y5 ta 1 j ds, y > C 


Example 4.3. Let Y be a rv defined in Example 14.11 Consider W = 4 Y and Z 
Using (14.41) the pdf of rv 1U is 


fw(w) = ^I r n 


-771+1,0 

o- — 

2 2 w 

ifli “1” bi 1? 1? 1)771+1,77 

m,n— m+1 

4nr=i^ 

^0, 25 {fli 1? 1? 1)771+1,77 


IL 
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(4.4) 

(4.5) 

(4.6) 


(4.7) 
so that 

□ 

= Y- 1 . 

'> 0 , 










































whprp k — TT m 1 FT'* r(ai+bi) 
where k lU=i ^ ai lU=m+i r( ai ) ■ 

Similarly by using (14.6[i the pdf of rv Z is 


fz(z) = 


(Ti 


UV2 

i= 1 




( 1) l)m+l,: 


( hi hi, 1, l)m+l,/i 


, z> 0 , 


Where k* = nZi)/i^nU,i r(oi) ■ 

Example 4.4. Let X be an exponential rv with pdf given by (14.51) . Consider Y = 6X and 
Z = X% where 6 > 0. Using (14.4ft the pdf of rv Y is 


y 


P 

O 

i_ 


Mv) = i* 0 

where /J = A 0. Now by using (14.51) the pdf of rv Z is 

fz(z) = 4 -/ 1 ? 

A® 


, y > o, 




_1 



z > o. 


Thus, Y is exponentially distributed with parameter /?, whereas Z follows the Weibull 
distribution (|4.6[) with parameters A, 0 (see |l2|). 

Theorem 4.3. Let X \and X 2 be two independent /-function rv’s with pdf’s given by (14.ip 

is given b? 

hi) A 


for l — 1,2. Then, the pdf of the rv Y = is given by 


fv(y) = ki% 


5y 




-1 -1 

^ a 2 (J _ 

where k = A_i 2 _i , m = mi + n 2 , n = n 1 +m 2 , p = pi + q 2 , q = q 1 + p 2} 5 = 5^ 5 2 ' 

XlWl )X2\ (J 2 ) 

and the sequence of parameters are 

(^*5 ^)l,p ^li)l,ni (1 ^2 j 2(7 2 ^2j ? &2 -®2j 1 )l ,<72 (^li? 0*1 ^li)ni+l,pij 


(J)j, Bj 1 ftj') iq (6y,(7i ft lj ) 1,771]. (f ^2/ 2(72 -A22? ^2 -^-2ii ^ 2 i)l,p 2 (^1? 5 0*1 5 /^l? )mi+l 

Proof. From (13.61) , we have for y > 0, 

Ms/) = .M- 1 -/>)] 


91- 


tf 1 <5 2 


(*C+ZOO 


AhOi 1 s)X2(2a 2 1 - cr 2 x s) (Wi 1 5 2 ^ ) ds, 


Xi(cri l )X2{cr 2 1 ) 2-Ki J c _ it 

which on using the definition of the /-function (11.11) completes the proof. □ 

Example 4.5. Let Xi,X 2 be two independent half-normal rv’s (14. 7 j) with pdf’s 

1 


fxX x i) = 


\iV2n 


/h° 

J o,i 


1 

e*. 

1 

51 

_ 1 



, Xi > 0 , 


where the parameters Aj > 0 for i — 1,2. Using Theorem 14.31 the rv Y = X\jX 2 follows 
Cauchy distribution (j4.8[) with pdf 


M.) = 


0,5,1 

0,|,1 


, y > °, 
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where the parameter A 


Ai/A 2 (see 0). 


4.3. The distribution based on the product of two /-functions. Let X be a rv with 
pdf 


fx (x) 


t, Tm,n 

hj l 1 p,q 


5x° 


(<3j, ylj, 

(bj, Bj, f3j) i )(J 


rk,l 

-^U^V 


TjX h 


(q, C{ , ^i) 1 } n 

(A, Dj, Pj)i, v 


x > 0 , 


0 


otherwise, 


(4.8) 


where a > 0, p > 0 and 6,rj ^ 0. Also, (a u A h af) ltP , (b 3 , B j: fy) hq , (c h Q, 7i)i,„, (dj, Dj, pf) 1)V 
m,n,k,l must conform to those restrictions in the definition of the /-function (11.11) . such 
that fx (x) > 0. Using Lemma [3. 11 the normalizing constant k\ is 


? 


k 


— 1 ® jn-\-k,m-\-l 

V 

(cijCj^i)!,/ (l-Zj -o- 1 B jl pa 1 B j ,{3 j ) 1>q (cj,Q, 7 i) /+ i jtt 

1 ^o -- 1 1 Q+U,P+V 

3. 

1 

b 

_ 1 

(dj , Dj ? Pj ) 1 ? /c (1 cii ^ A-i, ga)i,p (dj 5 Dj ? p^)fc+i,-u 


Special cases, (i) The generalized gamma distribution associated with Bessel function (see 
) has density 


fx(x) = 


T(a)ea 

-- 7 - 1,0 
— CLG a jlq 1 


x a e ax 0 F 1 (—; a; Xx ), x > 0, 


ax 


(a-1,1,1) 


7 - 1,0 

J 0,2 


— Xx 


( 0 , 1 , 1 ) (1 — a, 1 , 1 ) 


(4.9) 


where a, a > 0 and A ^ 0. 

(ii) The non-central chi-square distribution (put a — |,a = \ and A = | in (14. 9 j) ) is 


fx(x) = 


1 -—i 


2 iT( 


v f3x 


x 2 e 2 0 F 1 \ , x > 0, 


e 2 7 - 1,0 
o J o,i 


X 


7 - 1,0 

(3x 


2 

(l- 1 , 1 , 1 ) 

J 0,2 

- T 

(0,1,1) ( 1 -1,1,1) 


, (4.10) 


where v G N is the degree of freedom and j 3 > 0. 

(iii) The generalized gamma distribution associated with confluent hypergeometric function 
(see Q) has density 

fxix) = a ^- Y - x a ^e~(a+i)x iFl (£ ; Qj; x ), X > 0 , 

r («) 


a&(a + l) 1 ^ 1 


r(/3) 

where a,/3,a> 0. 


7-i,U 

J o,i 


(a + l)x 


(a-1,1,1) 


A; 


1,1 


— X 


( 1 -/ 3 , 1 , 1 ) 

(0,1,1) (1 — a, 1,1) 

(4.11) 


Remark 4.1. For k = u = v = 1, l = p = 0, \rj\ <1 and (ci, CV, 71 ) = (1,1,1), (di, £> 1 , pi) = 
(0,1,1) the pdf (14.81) reduces to (14. ip . 
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5. The /-Function inverse Gaussian Distribution 


Consider a rv X > 0 with pdf 
fxix) = < 


j^^^a—l^—ax—bx 1 jm,n 


Sx° 


( a u Ai, Q!i)i,p ^ n 

’ X>U ’ 


0 


with the following restriction on parameters: 


otherwise, 


(5.1) 


a > 0, b > 0 if a > 0, 

a > 0, b > 0 if a < 0, 

a > 0, b > 0 if a = 0. 

Further, a > 0,<5 ^ 0 and (a,, Ai, cu)i, p , (bj, Bj, /3j)i tq ,m,n must conform to those restric¬ 
tions in the definition of the /-function (11.11) . such that fx(x) > 0. In order to obtain the 
normalizing constant /c 2 of the /-function inverse Gaussian (I-FIG) distribution, we first 
determine its Mellin transform. 


Theorem 5.1. The Mellin transform of the I-FIG distribution is 

M fx (s) 




1=0 


J — Q — S + l J-m+1,71+1 

' 5 

(2 — a — s, cr, 1) (a*, A, <* 1 ) 1 , P (2 + Z - a - s, a, 1) 

S a 1 p+2,q+l 

a a 

(2 — a — s, a, 1) {bj^B^Pj) 


_ lck+s— 1 j-m+l,n+l 
U 1 p+l,q+2 


b a S 


(2 OL S, <7, 1) (Oj, Ai, Olf)ip 

(2 - a - s, a, 1) (1 - l - a - s, a,l) 


(< ab) l 

~i r 


Proof. Using the definition of Mellin transform and assuming the conditions for absolute 
convergence of integral involved, we have 


M fx {s) = h X a+s-2 e -ax-bx fjrn,n 


= ko 


ot+s—2—ax—bx 1 
X C- 


'0 


b 

H 

i_ 

( [C^i 5 ^ 4-25 ^z) l,p 

-Bj, 




dx 


- / x( w )(dx a ) w dw > dx 

2iri 


h_ I 

Jc u o 


a+s—aw—2—ax—bx 
X c- 


c 

-i 


dx > w dw. 


(5.2) 


Now 


^a.-\-s—crw—2—ax—bx 1 r> 

x e ax = 2 


ol-\-s—<jw —1 


Zfa+s-o™-i (2\/a6), (5.3) 


where K v (x) is the modified Bessel function of the second kind, also called Macdonald 
function, defined by 

= ,(/ (*)-/,(*)) (5 . 4) 
2 sm zz7r 

where G(t) is the modified Bessel function of the first kind with the following series rep¬ 
resentation: 


Ur) = ]u 


^ r(; + i/ + i)i! V 2 / 




(5.5) 
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Using fl5.3p - 05.5p and gamma reflection principle in (15.21) . we have 

k'2 


M fx (s) = 


2ni 


'c 


E 

1=0 


a— s+1 


r(a + s — aw — 1)T(2 — a — s + aw)x(ui) ( S 

a a 


T(l — o — s + aw + 2) 


b a+s-i r(a + s-aw- l)r(2 - a - s + aw)x(w) 

T(l + a + s- aw) 

Finally, by using the definition of the /-function (II.ip . the proof is complete. 
Corollary 5.1. The r-th moment for the I-FIG distribution is 


( ab) l 


dw. 


□ 


dr 


k 2 J2 

1=0 


) -ct-r Tm+l,n+l 

1 

On 

(1 - 0 - r, a, 1) (a*, A h Oi)i, p (1 + l - a - r, a, 1) 

S u 1 p+2,q+l 

a a 

(1 — a — r, a, 1) (bj, Bj, f3j) hq 


ra+r rm+l,n+l 
U 1 p+l,q+2 


b a S 


(1 H 1" 1 1) ^i)l,p 

(1 — o — r, < 7 ,1) (bj, Bj, /3j) ltq (-1 - a - r, a, 1) 


(ab) 1 


Hence, the normalizing constant of the I-FIG distribution is given by 


k 2 1 = 


1=0 


) —a jm+l,n+l 

' 5 

(l CM.) (J j 1) (^Z5 A-i) Q^z)l ? p (1 H - l CM) &) 1) 

S a 1 p+2,q+l 

a a 

(l-a, a, 1) (bj,Bj,fdj) i i9 


r- 



ra rm+l,n+l 
U 1 p+l,q+2 


b a 5 


(1 O, (7, 1) (<2j, 

(1 — o, a, 1) (bj.Bj^j)^ (-l-a,a, 1) 


(ab) 1 


Corollary 5.2. The Laplace transform of the I-FIG distribution is 

OO 

4W = 


J 1 j-m+l,n+l 

h 

(l-o, a,l) 


| (a + r) a p+2 ’ ,?+1 

(a + r) cr 

(1 -0,(7, 1) 

(6j, -Bj, f3j)l,q 


1=0 


1 a Tm+l,n+l 
U 1 p+l,q+2 


b a 5 


(1 O, (X, 1) (<2j, Oti)l,p 

(1-a, (7,1) (bj,Bj,f3j) hq (—1 — a, a, 1) 


((a + r)b) 
1! 


Special cases. The generalized inverse Gaussian distribution studied by Barndorff-Nielsen 

s. 


fx(x) = 




2 7 2/i Q (v^7) 

ft 

2 7 f K a (y/Fy) 


x a-l e ~\{Px+^x J ) 





(1.1.1) 

(J 

(0,1,1) 


x > 0, 


(5.6) 


where |h| < 1 ,/3, 7 > 0,o E M. Therefore the inverse Gaussian, inverse gamma and 
hyperbolic distributions follow as particular cases. 

(i) The inverse Gaussian distribution ^put a = — |,/3 = -=2 and 7 = A in (15.61) j. 


fx(x) = 


X AQ-m ) 2 


2ttx 3 


2/x z x 



A A _i_i_4 t 

= \l — e^x 2 e ^ 
2 t r 


T*” 1 7-1,0 

A 

( 1 , 1 , 1 )' 

J l,l 

U 

(0,1,1) 


x > 0, 
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where |h| < 1 and A,/i > 0. Note that K_i(x) = y|e X x 2 . 
(ii) The inverse Gamma distribution. 

X e 


fx(x) = 


m 

\ e 

m 




x~ e ~ 1 e~ Xx 17-1,0 


ri,u 

-' 1,1 


(1,1,1) 

(0,1,1) 


x > 0, 


where |6| < 1 and 6, A > 0. 

Remark 5.1. When 6 = 0, the I-FIG distribution in (15. lh belongs to the class of distribu¬ 
tions defined by the product of two /-functions (14. 8 j) as follows. 


fx{x ) = 


ax 



rm,n 

Sx a 

? -4. ii &i)l,p 

(a-1,1,1) 

1 P,Q 

(bj) Bj , Pj)l,q 


, x > 0, 


where 


U~ 1 _ -ajm,n+l 
^* ~ U I p+l,q 


(1 Cl, (T, 1) (dj, -T 7 , £^i)l,p 

(bj, Bj, fij)l t q 


Therefore the distributions with pdf fl4.9jl - fl4.llji are also the particular cases of the I-FIG 
distribution. 
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